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Abstract 

A search for stop production in i?-parity-violating super symmetry has been per- 
formed in e^p interactions with the ZEUS detector at HERA, using an integrated 
luminosity of 65 pb~^. At HERA, the i?-parity- violating coupling A' allows res- 
onant squark production, e^d q. Since the lowest-mass squark state in most 
supersymmetry models is the light stop, t, this search concentrated on produc- 
tion of t, followed either by a direct i?-parity-violating decay, or by the gauge 
decay to bxf ■ No evidence for stop production was found and limits were set 
on X[^i as a function of the stop mass in the framework of the Minimal Super- 
symmetric Standard Model. The results have also been interpreted in terms of 
constraints on the parameters of the minimal Supergravity model. 
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1 Introduction 



Many extensions of the Standard Model (SM) require a new fundamental symmetry be- 
tween bosons and fermions, known as supersymmetry (SUSY) [1]. This symmetry, hy- 
pothesizing the existence of supersymmetric partners of the SM particles, with similar 
properties but with spin changed by one half, controls the divergent higher-order loop 
corrections to the Higgs-boson mass. Despite numerous searches for such new parti- 
cles, no evidence has been observed, indicating that supersymmetry, if it exists, becomes 
manifest at scales beyond the present experimental limit. Since SUSY involves so many 
parameters, different experimental techniques complement each other. 

One important quantum number in supersymmetry models is i?-parity (Rp)- Its conser- 
vation ensures the conservation of both lepton and baryon number. Most of the searches 
performed at colliders drew conclusions under the assumption of Rp conservation. Never- 
theless, the most general supersymmetric extension of the SM Lagrangian contains terms 
which violate Rp and some of the possible i?p-violating (fip) scenarios are compatible with 
the present experimental constraints. 

One of the most interesting consequences of ftp scenarios is the possibility of producing 
single SUSY particles (sparticles) at colliders. Electron-proton collisions at HERA are 
well suited to the search for squarks, the scalar supersymmetric partners of quarks, since 
such states can be produced by an appropriate coupling of the incoming lepton and a 
quark in the proton. 

In most of the SUSY scenarios, the squarks of the third generation are the lightest; the 
present analysis is aimed of searching for the stop, t, the supersymmetric partner of the top 
quark. At HERA, the stop can be produced resonantly via e~^d — t, up to the ep centre- 
of-mass energy ^/s ^ 320 GeV. The stop decay can lead to distinctive topologies with 
a high-energy positron or neutrino and hadronic jets, which can be efficiently separated 
from the SM background. 

Direct searches for stop production have been already performed at HERA [2] , LEP [3-5] 
and Tevatron [6]. 

2 Stop phenomenology 

The i?-parity, defined as Rp = (_i)3-B+l+25^ jg multiplicative quantum number which is 
1 for particles and —1 for sparticles {B, L and S denote baryon number, lepton number 
and spin, respectively). i?-parity conservation would imply that supersymmetric particles 
are always pair produced and the lightest supersymmetric particle (LSP) is stable, a good 
candidate for cold dark matter. 
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If -R-parity is violated, it is possible to create single sparticles that decay to SM parti- 
cles [7]. The f^p terms in the SUSY superpotential are given by: 

W^^ = X^,kUL,^k + K,kL^QPk + K,kU^pk, (1) 

where the subscripts i,j,k are the generation indices, Li denotes the SU{2) doublet lepton 
superfield, Ei the SU{2) singlet lepton superfield, Qi the SU{2) doublet quark superfield 
and Di and Ui the SU{2) singlet down- and up-type quark superfields. The dimensionless 
Yukawa couplings Ajj^, X'^f., X'-jf, are free parameters of the model. 

In the case of stop production in ep collisions, the only terms involved are those param- 
eterized by the Yukawa coupling X'^^^. The partners of the left- and right-handed top, 
and tR, can mix together in two mass eigenstates, ti and t2, which, because of the large 
top mass, are usually strongly non-degenerate. Due to the chiral properties of the SUSY 
superpotential (Eq. 1) only the state contributes to the stop production cross section. 
In the regime of narrow- width approximation (NWA), the production cross section for 
the lighter state ii is hence: 

a{e+p^k) = ^ {coserX[,,r d{Ml/s,Ml), (2) 

where d{x, Q^j^ is the parton density of the d quark in the proton and 6^ is the mixing 
angle between ii and ^2- In this paper only the lighter stop, ii, denoted hereafter as 
t, has been considered since contributions from t2 are negligible for all the considered 
scenarios. The effects of the initial state photon radiation decreases the stop production 
cross section of ~ 5% (20%) for a stop mass of 150 (280) GeV and have been taken into 
account using the Weiszacker- Williams approach [8]. The NLO QCD corrections have 
been also included [9], they increase the LO cross section of ~ 20 — 25%. 

The t decays considered in this study are the -Rp-violating channel t —>■ e^d and the Rp- 
conserving decay i ^ Xi^i where Xi is the lightest chargincQ. These two channels provide 
a sufficiently large total branching ratio over all of the considered SUSY parameter space. 
The channel i — > y^t, where Xi is the lightest neutralinc@, contributes only at the highest 
stop masses, and even then it is below 10%. Branching ratios involving a heavier chargino 
or neutralino are small in most of the considered parameter space. 

* The variables x, and y, which is used later on in the paper, are the three Lorentz-invariant quantities 
characterizing the DIS processes. is the negative of the four-momentum-transfer squared, x the 
fraction of the proton momentum carried by the struck quark and y the inelasticity. 
The superpartners of the charged SU (2) gauge bosons and of the charged Higgs bosons mix together 
in two mass eigenstates named charginos. 

The superpartners of the neutral SU{2) gauge bosons and of the neutral Higgs bosons mix together 
in four mass eigenstates named ncutralinos. 
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The considered decays, including the cascade from the xti ^ire illustrated in Fig. 1. In 
the present paper only the hadronic W decays are considered, hence final states involving 
one positron with one jet (e-J), or more than one jets (e-MJ), and one neutrino and 
multiple-jets (z/-MJ) are studied. 

The results have been interpreted in the context of two different SUSY scenarios: the 
Minimal Supersymmetric Standard Model (MSSM) and the minimal Supergravity model 
(mSUGRA). 

In the unconstrained flp MSSM [10], the branching ratios for stop decay as well as the 
masses of the neutralinos, charginos and gluino^ are determined by the following MSSM 
parameters: the mass term /i, which mixes the Higgs superfields; the soft SUSY-breaking 
parameters Mi, M2 and M3 for the U{1), SU{2) and SU{3) gauginos, respectively; tan/3, 
the ratio of the vacuum expectation values of the two neutral scalar Higgs fields; and 
the ftp Yukawa couplings. The search for the stop was performed in the mass range 
100 — 280 GeV. In order to reduce the number of free parameters, the assumptions listed 
below were made: 

• no mixing between ti and {9i = 0); 

• only the Yukawa coupling was assumed to be non-zero; 

• SUSY scenarios in which the lightest neutralino Xi is not the LSP or is lighter than 
30 GeV, already excluded by LEP result [11], were not considered; 

• the GUT relations [12] for the gaugino mass terms Mi, M2 and M3: 

5 

Ml = -tan^Ow ■ M2, 
3 

M3 = ^sin^^vy-Ms 

C^EM 

were assumed. As a consequence, the gluino is always heavier than t, so that the decay 
t ^ t g is kinematically forbidden; 

• all the other sfermions (apart from the lighter stop) were assumed to have large masses 
that were fixed at 1 TeV. 

The total branching ratio due to the considered decay channels is > 70% for the range of 
parameters < 300 GeV and M2 between 100 and 300 GeV that was considered in this 
analysis. 

Figure 2 shows the branching ratios for four representative points in the H-M2 plane 
involving very different masses for and Xi ^ind for A'131 values equal to the limit given 
in Section 6.2. The shaded band indicates the sum of the three considered channels for 

^ Supersymmetric partners of gluons. 
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2 < tan/3 < 50. The total branching ratio and the contribution of the three channels are 
shown for the value tan (3 = 6. When the xt niass is larger than the stop mass, the fip 
decay (e-J) is the dominant channel; in the other cases the z/-MJ channel is generally the 
most relevant. 

In the mSUGRA [13], the number of free parameters is further reduced by assuming 
two universal mass parameters at the GUT scale, mo and mi/2, for all the sfermions 
and for all the gauginos, respectively. Radiative corrections are assumed to drive the 
electroweak symmetry breaking (REWSB), leading to consistency relations that allow the 
complete model to be fixed, based only on mo, mi/2, the sign of /i, tan (3, and the common 
trilinear coupling Aq. In the range of parameters used in this analysis, mi/2 < 180 GeV 
and mo < 300 GeV, the total branching ratio of the considered channels is in the range 
0.4 — 0.8 for mo < 200 GeV. For larger values of mo it decreases rapidly, since decay 
channels involving the heavier chargino, which are not considered in this study, become 
important. 



3 Data sample and experimental set-up 

The data used in this analysis were collected in the years 1999-2000. The total integrated 
luminosity was 65.1 ± 1.5 pb~^ of e'^p collisions. The proton and positron energies were 
^beam _ QgY ^nd E^^^^^ = 27.5 GeV, respectively, leading to a centre-of-mass energy 
of 318 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [14]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [15], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the 
polar-angl^ region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks is a{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [16] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 



^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18 /\/E for electrons 
and a{E)/E = 0.35//E for hadrons, with E in GeV. 

The luminosity was measured using the bremsstrahlung process ep —>■ epy. The resulting 
small-angle energetic photons were measured by the luminosity monitor [17], a lead- 
scintillator calorimeter placed in the HERA tunnel a.t Z = —107 m. 

4 Monte Carlo simulation 

The signal processes were simulated with the Monte Carlo (MC) event generator SUSY- 
GEN 3 [18]. It uses the exact matrix element for the production and for the decays 
of sparticles, includes initial- and final-state radiative corrections and is interfaced to 
Pythia 6.2 [19] for the hadronization of the final state. The program Suspect 2.1 [20] 
was used to solve the REWSB consistency relations that determine the sparticle mass 
spectrum at the electroweak scale in the mSUGRA model. 

The dominant background process to the e-J and e-MJ channels is neutral current deep 
inelastic scattering (NC DIS). For the e-J channel, the background is due to 2 ^ 2 scatters 
between high-x quarks and the positron. Backgrounds to the e-MJ channel occur in NC 
DIS events where multi-jet final states result from higher-order QCD effects. 

The primary background to the i^-MJ channel comes from charged current deep inelastic 
scattering (CC DIS) with multiple jets from QCD radiation. An additional background 
source involves photoproduction events for which the measured transverse momentum is 
large due to energy mismeasurement. 

The NC and CC events were simulated using the Heracles 4.6.1 [21] program with the 
Djangoh 1.1 [22] interface to the hadronization program and using the CTEQ5D [23] 
set of parton distribution functions (PDFs). In Heracles, corrections for initial-state 
and final-state electroweak radiation, vertex and propagator corrections, and two-boson 
exchange are included. The colour-dipole model of Ariadne 4.10 [24] was used to simulate 
the order as plus leading-logarithmic corrections to the quark-parton model. The MEPS 
model of Lepto 6.5 [25] was used as a systematic check. Both programs use the Lund 
string model of Jetset 7.4 [26] for the hadronization. 

Photoproduction events were simulated using the Herwig 6.100 [27] generator, using the 
CTEQ4L [28] proton PDFs. Both direct and resolved photoproduction were considered. 
In the direct case, all of the photon energy participates in the hard scattering, whereas, 
for the resolved process, only a fraction of the photon energy, associated with a parton 
constituent of the photon, participates in the hard subprocess. For the simulation of the 
resolved subprocess, the GRV-G [29] photon PDFs were used. 
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5 Event selection 



The signal events are characterized by a high-energy lepton in the final state. In the case 
of a positron in the final state (e-J and e-MJ channels), the trigger selection was based 
on a standard neutral current trigger, which required a scattered positron, as used in 
searches for resonance states decaying to eq [30] and in ZEUS NC DIS studies [31]. For 
the neutrino case (z/-MJ channel), a trigger selection based on a missing- requirement 
and already employed in the ZEUS CC DIS analysis [32] was used. 

The offline signal-search procedure was performed in two steps [33]. Initially, a preselection 
was applied to select NC or CC events. Finally, more restrictive selections, optimized to 
get the best limits in the case of no signal, were applied. 

5.1 Preselection for e+ final states 

The following conditions, some of which were also used at the trigger level with a lower 
threshold, were designed to select a sample of high-Q^ NC events: 

• Z-coordinate of the event vertex compatible with an ep interaction, \Zyt^\ < 50 cm; 

• a high-energy positron reconstructed from calorimeter and tracking information [34]. 
A positron energy Eg > 8 GeV was required. This cut was increased to Pt^e > 20 GeV 
(Pr,e is the transverse momentum of the positron measured by the calorimeter) for 
very forward positrons {6^ < 0.3, where 6e is the positron polar angle) which are 
outside the acceptance of the CTD; 

• 45 < E — Pz < 70 GeV, where E and Pz are the total energy and the Z-component 
of the total momentum of the final state. For NC DIS events, where only particles in 
the very forward direction escape detection, E — Pz ^ 2E^'^^^ = 55 GeV; 

• Qda > 1000 GeV^ and 0.2 < ^da < 0.98, where Qda ^^'^ 2/da are the DIS kinematic 
variables reconstructed using the double angle method [35]. The above conditions 
were imposed in order to restrict the search to a region where the signal is enhanced 
with respect to NC DIS and the reconstruction of the kinematic variables is reliable; 

• Mex > 100 GeV, where Mex is the invariant mass of the positron and the hadronic sys- 
tem evaluated using the following relation that exploits the conservation of momentum 
and E - Pz: 



The sum runs over the final-state positron and all other energy deposits with a polar 
angle > 0.1, to exclude contributions from the proton remnant. 




(3) 
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The bias and resolution of the reconstructed mass were evaluated using the signal MC. 
On average, the mass was slightly overestimated at low masses (3% at 100 GeV), while 
the agreement improved towards high masses (< 1% above 150 GeV). The resolution 
varied between 5% and 1.5% in the mass range 100 — 280 GeV. After the preselection 
cuts, 2368 events remained, in good agreement with the expectation of the SM MC of 
2430I252 events, where the error is dominated by the systematic uncertainties described 
in Section 16.11 The SM prediction is dominated by the NC DIS contribution. 

Figure 3 shows the distributions of Pt.c, Vda, IoSioI'^da) -PT,antipar/-PT,par for data and 
MC; reasonable agreement is seen for all variables. The quantities Pr.par and Py.antipar are 
the parallel and antiparallel components of the hadronic transverse momentum (Pr.had) 
defined as: 



where the sums are over calorimeter deposits with polar angle 6 > 0.1, excluding the 
identified positron. The ratio -PT,antipar/-Pr,par is used in the final selection to separate 
one-jet events (PT,antipar/-Pr,par ~ 0) from multi-jet events (PT,antipar/-Pr,par > 0). 

5.2 Preselection for u final state 

Events with a neutrino in the final state have a topology similar to CC DIS. The following 
selection cuts were applied in order to select a sample of high-Q^ CC events and suppress 
the non-ep contribution: 

• Z-coordinate of the event vertex compatible with an ep interaction, |Zvtx| < 50 cm; 

• no reconstructed positron satisfying the same criteria used in final-state preselec- 
tion; 

• high missing transverse momentum, PT,miss > 20 GeV, where Pr.miss is the missing 
transverse momentum as measured by the CAL; 

• 0.2 < yjB < 0.95, where yjB is reconstructed using the Jacquet-Blondel method [36]. 

The analogue of Eq. 3 for the invariant mass of the z/-hadronic system was derived as- 
suming that the missing Px and E — Pz resulted from a single neutrino: 




The mass resolution varied between 10% and 3% in the mass range 100 — 280 GeV. On 
average, the mass was shghtly underestimated at high masses (1.5% at 280 GeV), while the 
agreement improved towards low masses (< 1% above 120 GeV). After the CC preselection 
cuts 265 events survived, in good agreement with the expectation of the SM MC of 
The SM prediction is dominated by CC DIS events, with a small contribution coming from 
photopro duct ion processes. 

Figure 4 shows the distributions of Py^miss, Ujb, IoSioI'^jb)' where the Qj^ is recon- 
structed using the Jacquet-Blondel method, and -Pr,antipar/-Pr,par for data and MC; rea- 
sonable agreement is observed for all the variables. 



5.3 Final selection for e+ final state 

The final selection for the channels with a final-state positron was designed to reduce 
further the contamination from NC DIS by requiring high-Q^ and high-y events. The 
following cuts were applied: 

• > 3000 GeV^; 

• Uda > Ucut, where i/cut was optimized as a function of the reconstructed mass using the 
SM MC and ranges between 0.7 and 0.4 for masses between 100 and 280 GeV. This 
cut exploits the different y-dependence of NC DIS, steeply decreasing as y"^, and of 
a scalar resonance, which has a substantial contribution from large y. 

Finally, a cut on -Pr,antipar/-Pr,par was used to produce two samples enriched with either 
one-jet or multi-jet events. The e-J (e-MJ) final sample was defined requiring: 

• -Pt,antipar/-Pt,par < (>) 0.05. 

Signal efficiencies were evaluated by generating samples of signal events using Susygen 
for different values of the MSSM or mSUGRA parameters. For the e-J channel, the 
efficiencies ranged between 10 and 45% in the mass range 100 — 260 GeV, decreasing to 
20% at 280 GeV. For the MSSM scenario, the efficiencies for the e-MJ channel were in the 
range 5 — 25% for stop masses between 200 and 280 GeV, depending mainly on the masses 
of stop and xt'i efficiency decreased towards lower and higher stop masses. For the 
mSUGRA scenario, the efficiencies for the e-MJ channel were in the range 5 — 15% in most 
of the parameter space. Table 1 shows good agreement between the number of selected 
events and SM expectation. Figures 5a and 5b show reconstructed mass distributions 
of data and SM MC for the e-J and e-MJ preselection and final samples. The data 
distributions are well described by the SM simulation. 



8 



5.4 Final selection for v final state 

In order to enhance stop sensitivity and reduce further the contribution of CC DIS, the 
final selection required: 

• yjB > 0.6; 

* Pt, 'Aatipax / Pt, pax ^ 0.1. 

For the MSSM scenario, the efficiencies were in the range 15 — 35% for stop masses 
between 180 and 280 GeV, depending mainly on the masses of the stop and Xii the 
efficiency decreased towards lower and higher stop masses. For the mSUGRA scenario, the 
efficiencies were in the range 5 — 20% in most of the considered parameter space. Table 1 
shows good agreement between the number of selected events and the SM expectation. 

Figure 5c shows reconstructed mass distributions of data and SM MC for the i^-MJ pre- 
selection and final samples. The MC simulation also in this case describes the data 
reasonably well. 



6 Results 

Since no evidence for stop production was found, limits at 95% CL were set using a 
Bayesian approach. The limits were set for two different SUSY scenarios: the uncon- 
strained MSSM model and the mSUGRA model (see Section 2). 

6.1 Systematic uncertainties 

In the calculation of the upper limit on A'^g^, several sources of systematic uncertainties 
were considered. The following systematic uncertainties on the SM background expecta- 
tion were evaluated: 

• the uncertainty from the proton PDFs, evaluated using the procedure suggested by 
the CTEQ group [37], was ±4% for z/-MJ and ±2% for e-MJ and e-J; 

• the uncertainty on the calorimeter energy scale of ±1% (±2%) for the electromagnetic 
(hadronic) section led to an uncertainty on the SM event rate of ±5% for i^-MJ and 
+3^° for e-MJ and e-J; 

• the use of MEPS instead of Ariadne to simulate the QCD cascade led to an uncer- 
tainty of —3% for I/-MJ and —6% for e-MJ and e-J; 

• the uncertainty in the integrated luminosity measurement was ±2.25%. 

In addition, the following uncertainties related to the signal simulation were considered: 
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• the uncertainties in the signal efficiency due to interpolation between different SUSY 
scenarios was ±15%; 

• the theoretical uncertainty on the signal cross section due to the uncertainty in the d- 
quark parton density [37] in the proton varied from ±3% to ±80% for masses between 
100 and 280 GeV. 

6.2 Limits for the MSSM model 

Assuming the MSSM model, the upper limits on A'^g^ were evaluated as a function of the 
stop mass. A scan of the mass spectrum in 1 GeV steps was performed using a sliding 
window of ±2(Jmj for Mix < 250 GeV (£ = e or z/), where (Tmj is the stop mass resolution. 
For masses larger than 250 GeV, where the SM background is smaller and the expected 
signal width larger, the condition M^x > — 2aM^ was apphed. 

At each stop mass, the 95% CL limit on A'^g^ was evaluated using, for each channel, 
the data events, the SM predictions and the signal expectation for the corresponding 
mass window. The signal cross section was calculated in the NWA (Eq. 2), including 
initial-state radiation for the incoming positron [8] and the next-to-leading-order QCD [9] 
corrections, using the CTEQ6 [37] set of parton densities, while the branching ratios for 
the different channels and MSSM scenarios were taken from the Susygen simulation. The 
total likelihood was evaluated as the product of the Poissonian likelihoods of each channel. 
The systematic uncertainties described in Section 6.1 were included in the likelihood 
function assuming Gaussian probability densities. A Bayesian approach assuming a fiat 
prior for the signal cross section was then used to produce the limits. 

Figure 6 shows the 95% CL limit on \'^^^ as a function of the stop mass for the range 
-300 < /i < 300 GeV, 100 < Ms < 300 GeV and 2 < tan/5 < 50. The limits for masses up 
to 250 GeV improve on the low-energy constraints from atomic parity-violation (APV) [38] 
measurements (dashed line) and do not depend strongly on the different SUSY scenarios. 
The HI collaboration obtained similar constraints [2] using similar SUSY scenarios. 

6.3 Limits for the mSUGRA model 

For fixed values of X'l^i-, constraints on the mSUGRA parameters can be set in the plane 
(mo,mi/2), when tan/9, Aq and the sign of are fixed. The parameter Aq enters only 
marginally at the electroweak scale and was set to zero. Limits at 95% CL were evaluated 
using a scan of the reconstructed mass spectrum and the same Bayesian approach as in 
the MSSM case. 
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Figure 7 shows the 95% CL excluded area in the plane (mo,mi/2) for A'^3^=0.3, tan/3 = 6 
and /i < (hatched area). The dark region corresponds to values of parameters where no 
REWSB solution is possible, while the light region corresponds to neutralino masses less 
than 30 GeV, already excluded by LEP [11]. The dashed lines indicate curves of constant 
stop mass close to the border of the excluded area. Stop masses can be excluded up to 
250 GeV for mo smaller than 240 GeV. The effects of the SUSY radiative corrections on 
the sparticle mass spectrum is included in SUSPECT and have been taken into account. 
Such effects increase the stop mass and consequently worsen the limits especially at large 
rriQ. For example the point mo = 200 GeV, mi/2 = 110 GeV is at the boundary of the 
ZEUS exclusion region and corresponds to = 256 GeV. The same point, if SUSY 
radiative corrections are neglected, corresponds to = 243 GeV and would be well 
inside the ZEUS excluded region. 

A scan towards large tan/3 was performed assuming M = mo = mi/2- Figure 8 shows 
the limits on M as a function of tan/3 for X[^i=0.3 and fi < 0. The limit on M slightly 
increases from 130 to 140 GeV in the range 6 < tan/3 < 40. For larger values, it drops 
because the large mixing in the f sector results in a light fi state into which the t can 
decay. The efficiency for detecting such decay is low. The effect of the SUSY radiative 
corrections is to slightly decrease the overall limit and to shift towards larger tan/3 the 
point where the stau branching ratio opens up; neglecting radiative corrections the limit 
drops at tan /3 ~ 37. The HI collaboration obtained comparable constraints [2] using the 
same mSUGRA scenarios. 

6.4 Comparison to results from other colliders 

Studies on stop in flp SUSY scenarios have been performed both at LEP [3-5] and the 
Tevatron [6], looking for the production of stop pairs. LEP mass limits for the stop, in 
the case of A' > 0, were obtained by the OPAL [3] and ALEPH [4] collaborations and are 
in the range 85 — 98 GeV. The CDF collaboration [6] set a stop mass limit at 122 GeV 
assuming X'^.^f^ > and a branching ratio B{t hr) = 1. 

A more interesting comparison between HERA and Tevatron sensitivities can be done 
by looking at Tevatron results for leptoquark (LQ) production. The DO collaboration 
published limits on leptoquark masses as a function of the branching ratio B[LQ — »• 
eg) [39]. Since leptoquarks and squarks have analogous production mechanisms, such 
limits can be converted into limits on the stop mass as a function of X[^i [40, 41] and 
directly compared with the results of this analysis. In the case of the MSSM scenarios, 
DO limits are competitive with those of HERA only for the largest values of M2 and 
where the flp decay i ^ eq dominates due to the large chargino mass. For lower values 
of M2 or the gauge stop decays are relevant and the ZEUS limits improve over those 
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from DO for masses larger than 150 GeV. Figure 9 shows the comparison between ZEUS 
and DO hmits for three different regions of the unconstrained MSSM parameter space. In 
the mSUGRA scenarios considered here, the gauge stop decays are always relevant and 
thus the ZEUS limits are more stringent than those from DO. 

7 Conclusions 

A search for stop production in e^p collisions at HERA was performed using an integrated 
luminosity of 65 pb~^. No evidence was found for resonances in the decay channels with 
jet(s) and one high-P^ positron or neutrino. The results have been interpreted in the 
framework of the i?-parity-violating MSSM, setting limits on the Yukawa coupling \'^^^ 
as a function of the stop mass. These limits exhibit a weak dependence on the MSSM 
parameters /x, M2 and tan/3 and improve on limits from Tevatron in a large part of 
the considered parameter space, and on limits from low-energy atomic parity-violation 
measurements for stop masses lower than 250 GeV. Direct limits on the stop mass have 
also been derived within the mSUGRA model. In this model only five free parameters 
determine the full supersymmetric mass spectrum. In this case, for = 0.3, tan/3 = 6, 
/X < and = 0, stop with masses as high as 260 GeV are excluded for a large part of 
the parameter space. 
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Channel 


QIa (GeV2) 


2/cut(Mi) 


-Pr.antipar/ -Pf.par 


Data 


SM 


Eff. MSSM 


e-J 


> 3000 


0.4-0.7 


< 0.05 


85 


74.5^1^ 


0.3 


e-MJ 


> 3000 


0.4-0.7 


> 0.05 


63 


58.81^:° 


0.15-0.2 


I/-MJ 




0.6 


> 0.1 


19 


20.9t\i 


0.15-0.35 



Table 1: Summary of final selection cuts, number of observed and expected events 
for the SM and signal efficiencies for MSSM (Mi = 220 GeV) for the different 
channels discussed in the text. 
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Figure 2: Branching ratios as a function of the stop mass for different fi and 
M2 values, and for X[^i equal to the ZEUS limit for the considered scenarios. The 
three channels considered in the analysis and their sum are shown for tan/5 = 6. 
The shaded band is the range of the total branching ratio for 2 < tan/3 < 50. 
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Figure 3: Comparison between data (dots) and SM MC (histograms) for (a) Pt^e, 
(b) uda, (c) logio(Ql)A) (^iT'd (d) PT,anttpar/ PT,par after the preselection for e+ final 
state. The band represents the SM expectation with its uncertainty. The MSSM 
signal for the e-MJ channel (dashed line) with = 220 GeV , M2 = 100 GeV and 
fj, = —300 GeV is also shown, with arbitrary normalization. 
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Figure 4: Comparison between data (dots) and SM MC (histograms) for (a) 
Ptmss, (b) yjB, (c) logio(QjB) (d) Pt,antipar/ Pt,par after the preselection for 
V final state. The band represents the SM expectation with its uncertainty. The 
MSSM signal (dashed line) with = 220 GeV, Ms = 100 GeV and fi = -300 GeV 
is also shown, with arbitrary normalization. 
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Figure 5: Reconstructed mass for (a) e-J, (b) e-MJ and (c) v-MJ channels. The 
data (dots) and SM MC (histograms) after preselection (empty circles and light 
histograms) and final selection (filled circles and dark histograms) are shown. The 
hand represents the SM expectation with its uncertainty. 
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Figure 6: Exclusion limits on X[.^^ as a function of the stop mass for the MSSM 
model. The light (dark) region is excluded in all (part of) the considered SUSY 
parameter space. The region above the dashed line is excluded by low- energy atomic 
parity-violation (APV) measurements. 
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Figure 7: Exclusion limits for mSUGRA with tan(3=6, X[^i=0.3, /x < and 
Ao = (hatched area). The dark-shaded region corresponds to values of parameters 
where no radiative electroweak symmetry- breaking solution is possible. The light- 
shaded region corresponds to neutralino masses (LSP) less than 30 GeV , already 
excluded by LEP results. The dashed lines indicate the curve of constant stop mass 
close to the border of the excluded region. 
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Figure 8: Exclusion limit for mSUGRA on the mass parameter M (M = 
rriQ = mi/2) as a function of tan j3 for \'^^^=0.3, fi < and Aq = 0. The dark- 
shaded region corresponds to values of parameters where no radiative electroweak 
symmetry-breaking solution is possible. 
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Figure 9: Comparison between ZEUS, DO and atomic parity violating (APV) 
limits in the MSSM scenario for tan (3 = 6 and (a) low M2, (b) low \fi\ and (c) high 
M2 and The regions excluded by ZEUS and DO are shown by the dark-shaded 
area and by the area above the full line, respectively. The regions excluded in part 
of the parameter space are the light-shaded area (ZEUS) and the area between the 
full and the dashed line (DO). The region of exclusion for atomic parity-violation 
(APV) is above the dotted line. 
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